The Anjials of Applied Probability 
2011, Vol. 21, No. 1, 115-139 
DOI: 10.1214/10-AAP689 

@ Institute of Mathematical Statistics, 2011 
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We consider waves propagating in a randomly layered medium 
with long-range correlations. An example of such a medium is stud- 
ied in [19] and leads, in particular, to an asymptotic travel time de- 
scribed in terms of a fractional Brownian motion. Here we study the 
asymptotic transmitted pulse under very general assumptions on the 
long-range correlations. In the framework that we introduce in this 
paper, we prove in particular that the asymptotic time-shift can be 
described in terms of non-Gaussian and/or multifractal processes. 

1. Introduction. Wave propagation in random media has been exten- 
sively studied for many years from both theoretical and applied points of 
view. In particular, the study of the effective shape of an acoustic pulse prop- 
agating through a layered medium has attracted a lot of attention [1,5, 24]. 
Recently, applications to time reversal [10] have also been the subject of 
much interest. Currently, there is also a strong interest in problems related 
to noise and correlations [11]. In all these situations the statistical properties 
of the medium are important since they affect the statistical properties of 
the wave field. 

In [5] the authors consider an acoustic pulse propagating in a 
one-dimensional random medium with rapidly decaying correlations. They 
rigorously prove the classical result of O'Doherty and Anstey [20] which es- 
tablishes that the effective transmitted pulse is characterised by determin- 
istic spreading and a random time-shift. More precisely, the deterministic 
spreading is expressed as a convolution with a Gaussian density, and the 
random time-shift is described in terms of a Brownian motion. 

More recently, wave propagation and also homogenization in random me- 
dia with long-range correlations and /or defined in terms of fractional Brow- 
nian motion [2, 13, 19, 23] have been considered. In [19], we extend the 
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result of [5] to such a framework. Then, the asymptotic description of the 
transmitted pulse is dramatically different from what happens in a mixing 
case. Indeed, the pulse keeps its initial shape, and its random time-shift is 
now described in terms of a fractional Brownian motion whose Hurst in- 
dex depends on the decay rate of the correlation function of the random 
fluctuations. We considered in [19] a particular form of a random process 
describing the medium, such that it was roughly speaking close to a Gaussian 
process. Thus, it still remains to study more general cases under the long- 
range assumption. This is the aim of the present work. We establish that 
under general long-range assumptions on the medium, the effective pulse 
still keeps its initial shape as observed in [19], but the time-shift can be very 
different, non-Gaussian, for instance, depending on the form of the random 
fluctuations. Moreover, our general result allows us to deal with media with 
a decay rate of correlations that varies along the propagation direction. This 
leads to an effective time-shift, described in terms of a multifractional ran- 
dom process which is, roughly speaking, a fractional Brownian motion with 
a varying Hurst index that reflects the nonhomogeneity of the propagation 
medium. 

In Section 2 we introduce the problem and review the basic wave decom- 
position approach. Next, we establish the general technical result (Theorem 
3.1) in Section 3. We apply this general result to non-Gaussian media in 
Section 4, to multifractal Gaussian media in Section 5 and to multifractal 
non-Gaussian media in Section 6 where we prove the main result of the pa- 
per (Theorem 6.1). We present a numerical illustration in Section 7. Finally, 
Section 8 is devoted to the derivation of Theorem 3.1. 

2. Preliminaries. 

2.1. Wave decomposition. The governing equations are the nondimen- 
sionalized Euler equations giving conservation of moments and mass 



where t is the time, z is the depth into the medium, is the pressure and 
the particle velocity. The medium parameters are the density and the 
bulk- modulus (reciprocal of the compressibility). We assume that is 
a constant identically equal to one in our nondimensionalized setting, and 
is modeled as follows: 



(2.1) 




dp^ 





(2.3) 




for z€ [0,Z], 
for z G M - [0, Z] 
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where ^'^ is a centered random process. The number e > is a parameter that 
aU quantities depend on. As we wiU see below it is introduced to describe 
the scales of the problem. 

We introduce the right- and left-going waves 

(2.4) A'=f + u'' and B^ = u^-p''. 
The boundary conditions are of the form 

(2.5) A^{z = {),t) = f{t/e^) and B^z = Z,t)=0 

for a positive real number r > and a source function /. In order to deduce 
a description of the transmitted pulse, we open a window of size e'^ in the 
neighborhood of the travel time of the homogenized medium and define the 
processes 

(2.6) a''{z,s)=A'{z,z + e''s) and b'{z,s) = B''{z,-z + e'^s). 

Observe that the background or homogenized medium in our scaling has a 
constant speed of sound equal to unity and that the medium is matched so 
that in the frame introduced in (2.6) the pulse shape is constant if ^u^ = or 
if we consider the homogenized medium [10]. We introduce next the Fourier 
transforms and of and ¥ , respectively, 



oo roc 



o^(z,w)= / e"^'a^z,s)ds and b^z,u}) = / e''^'b'{z,s)ds, 

J —oo J —oo 

that satisfy 

(2.7) ^ = '-^i^'izW - e-^^'^'/'^b^), a'{0,u) = f{oj), 

(2.8) — = —u^z){e^''^'/'^a'-b'), b^Z,oj)=0, 
dz 2 

where we use the notation 

(2.9) iy' = ^. 

Following [5, 10] we express the previous system of equations in terms of the 
propagator P^i^z) which can be written as 



(2.10) P^(z) 
and that satisfies 



/5S(2) a%[z) 



(2.11) ^(.)=7^^(^^,.)p^(z), P^(^ = 0)=(J \ 
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with 



iijj ^ , , ( 1 —€ 



Defining next the transmission coefficient and the reflection coefficient 
^5 by 

(2.12) T^(^)= 1 and i?^(z) = ^^^^ 



we can write 



(2.13) a-(Z, = ^ e---r^(Z)/(a;) 
and 

(2.14) 6^(0, s) = e-''^R'jZ)fiu)div. 

Hence, we shall study the asymptotics of the propagator in order to 
characterize and as e goes to 0. 

2.2. A short-range medium. We recall now what happens in a mixing 
(or short-range) model when r = 1 and /U^(z) = v^z/e^). We assume that 
1/ = ^om where <I> is a bounded function and m is a centered Markov process 
with an invariant probability measure whose generator satisfies the Fredholm 
alternative. This implies that the covariance function z i— )■ E[z^(0)z^(2:)] is 
integrable and then the correlation length a of the medium is finite 



poo 

= / E[i/(0)i/(z)]dz G [0,oo). 

JO 



This property is the mixing property or the short-range property. It is well 
known [5, 10] that under these assumptions the propagator equations 
converge to a system of stochastic differential equations driven by indepen- 
dent Brownian motions from which we can deduce that a^[Z,s) — >a{Z,s) 
as e goes to with 

(2.15) a{Z,s) = {f*G){s-B), 

where G is a centered Gaussian density with variance a^Z/2 and B a Gaus- 
sian random variable that can be expressed in terms of a Brownian motion W 
as B = aW{Z)/y/2. Proving this result involves using the Diffusion Approx- 
imation Theorem [10] to get an asymptotic propagator from which we can 
deduce the expression of the limit a{Z, s). Notice that, whereas the variance 
of B depends in particular on the result does not depend qualitatively 
on $ in the sense that B remains Gaussian whatever <I> is. 
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2.3. A long-range medium. In [19] , the propagation in a long-range medium 
is investigated. The model considered is defined in terms of a fractional 
Brownian motion. More precisely, we assume that z^*" has the form 

jy'^(z)=e^~^i^(^-^^ forzG[0,Z], 

where k> and is a process that is expressed as z^(z) = for every 

z where: 

• $ is an odd C°°-function; 

• m is a Gaussian process, centered, stationary and has a correlation func- 
tion which has the following asymptotic property as z goes to oo: 

(2.16) rmiz)=nmiO)miz)]^c„^z-\ 7e(0,l). 

The property (2.16) implies that the covariance function of v is not inte- 
grable 

|r,y(2:)| dz = oo, 

which means that the correlation length is infinite. This is the so-called long- 
range property. We mention that a typical example of a process satisfying 

(2.16) can be constructed as 

(2.17) m{z) = WH{z + l)-WH{z), 

where Bh is a fractional Brownian motion (fBm in short) with Hurst pa- 
rameter H > 1/2. 

We assume r, k and 7 satisfy t — k = -f. In this case, we proved that 
a'^{Z, s) — )• a{Z, s) with 

(2.18) aiZ,s) = f{s-B), 

where B a Gaussian random variable. We can write B as B = (T//W//(Z) 
where Wh is a fractional Brownian motion with Hurst parameter H = (2 — 
7)/2 and an is a positive constant that depends on H and 



3. Medium assumptions and main technical result. The results presented 
above show that the asymptotic behavior of the pulse shape a^{Z, s) strongly 
depends on the statistical properties of z^. The pulse shape is affected under 
short-range assumptions whereas it does not change under the long-range 
assumptions described above. In Sections 4 and 5 we carry out the analysis 
of the particular long-range media that we consider in this paper. To facili- 
tate this analysis we establish in this section a theorem under the following 
general assumptions on i^^ = /-i^/e^. Let A > and define: 
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• Assumption Ai: As e goes to 0, the finite-dimensional distributions of the 
process {J^ i''^{z')dz'}z converge to those of a process V = {V{z)}z with 
finite second-order moments. 

• Assumption A2 (A) : There exist two symmetric, continuous and two- variable 
functions j:[0,Z]^ [7-,7+] C (0,1) and R:[0,Z]^ (0,oo) such that 
for every 6 > 0, there exists zs > sufficiently large such that for every 
zi, Z2 and e satisfying l^i — 2:2] > e'^^zs, 



E[i^%z,)u%Z2)]-R{z^,Z2)\zi-Z2r'-''''''H<6R{z,,Z2)\zi-Z2r'^'''''l 



• Assumption A3 (A). For every p> there exist Cp> and 7p G (0, 1) such 
that |E[i/^(zi)z/'^(z2)]| < Cp\zi — Z2\~'~^f for every e > and \z\ — Z2\ < e^p. 

Assumption Ai corresponds to the convergence of the travel-times. Assump- 
tions A2(A) and A3 (A) are long-range assumptions for nonstationary pro- 
cesses. They describe how the long-range property varies with the propaga- 
tion distance. In particular, these enable us to apply the next theorem to 
multifractal media (Sections 5 and 6), which are nonhomogeneous. 

Here we give the main technical result of this paper. This theorem is 
next used in Sections 4, 5 and 6 to establish the asymptotic pulse shape 
respectively in non-Gaussian and multifractal media. 

Theorem 3.1. Assume that there exists A > such that Ai, A2(A) and 
A3(A) are satisfied. Then, as e goes to 0, {a^(Z, s)}<j converges in distribution 
in the space of continuous functions endowed with the uniform topology to 
the random process {a(Z, s)}^ that can he written as 



Theorem 3.1 establishes that, under general long-range assumptions, if 
the travel-time converges then the asymptotic pulse keeps its initial shape 
but its time shift is described in terms of the asymptotic travel-time. As 
recalled in Section 2.3 this fact was observed in a particular case in [19]. 
In fact, the result of [19] follows from Theorem 3.1. Indeed, the model pre- 
sented in Section 2.3 satisfies Ai, A2(2) and A3 (2). In particular, the finite- 
dimensional distributions of {Jq ^^{z') dz'}z converge to those of the process 
{2aHWH{z)}z, so that the asymptotic pulse is of the form (2.18). 

Notice that the framework we study in this paper is in dramatic contrast 
with the mixing case where we observe a pulse spreading in addition to the 
time-shift. This is not so surprising if we remark that Theorem 3.1 does not 
apply to a process defined as in Section 2.2 by i/^(z) = e~^i'{z/e'^) where 
1/ is a mixing process. Indeed, if such a process satisfied assumption A2(A) 
for some A > 0, then we would have 



(3.1) 



a{Z,s) = f{s-'^V{Z)). 
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for some c* > 0, z* G [0,Z] and 7* G (0,1), which contradicts the mixing 
assumption that gives 

POO 

E[i/^(z)i/^(0)]dz < / E[z^(z)i/(0)](iz<oo. 



z* JO 

To conclude this section we present a heuristic description of the hnk 
between the mixing and the long-range cases. For every uj we define 

where for every z £ [0,Z] by 

vl{z)= f\'{z')dz', 
Jo 

V2A^)= [ J^'{z')cos(2uj^^ dz', 



vl^{z) = j\'{z')s\n{2u:^^ dz' . 



In both cases the three-dimensional process plays a crucial role. In the 
mixing case converges to the three-dimensional (nonstandard) Brownian 
motion (-Bi, i?2,i^, ^3,1^)- In the proof of the convergence 

a'{Z,s)^a{Z,s) = {f*G){s-B) 

one then observes that the Gaussian variable B can be written as i? = 
B\{Z)/2, and that the Gaussian density G derives from B2^u and B^.^^ 
[5, 10, 12]. In the long-range case, let us assume that converges to the 
three-dimensional process (F, 0,0). This fact was already observed in [18] 
for the fractional white noise. Now if we substitute (Si, ^2,^^;, .63,1^) with 
{y, 0, 0) in the expression of the limit a{Z, s) we obtain B = V{Z)/2, G = 5q 
(because in fact (5=1) and hence a{Z, s) = f{s — V{Z)/2). This is what we 
establish in this paper, in particular by proving the convergence of and 
the substitution mentioned just above. 

4. Non-Gaussian asymptotics. In this section we study the case where 
i/^ has the form 

u'^(z)=e''-^u(^-^^ forzG[0,Z], 

where k > and is a process that is assumed to have the form 

i/{z) = ^{■m{z)) 

for every z where: 
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• <I> is a continuous function such that ^{ao x •) has a Hermite index equal 
to KgW, where = E[m(0)2]. 

• 771 is a continuous Gaussian process, centered, stationary and has a corre- 
lation function Vm which has the following asymptotic property as z goes 
to oo: 

(4.1) rmiz) = E[m(0)m(z)] ~ c^z~\ 
where < 7 < 1/K. 

We denote the Kth Hermite coefficient of <I?((To x •) by 

J{K)=EmaoX)PK{X)], 

where X ~ AA(0, 1), and Pk is the i^Tth Hermite polynomial. Applying The- 
orem 3.1 we get the following result. 

Theorem 4.1. Assume that t — k = jK. Then, as e goes to 0, {a^{Z,s)}s 
converges in distribution in the space of continuous functions endowed with 
the uniform topology to the random process {a{Z,s)}s that can be written as 

(4.2) aiZ,s) = fis-^Wlf{Z)), 

where is the Kth Hermite process of index H = [2 — ^K)/2 G (1/2, 1) 
defined for every z by 

K/2 ,. K 

(4.3) w^{z) = ^ GH,K{z.xi,...,XK)T\B{dxk) 

with 

g^,x(.,.„...,x^)= ^,c7(F)^Ef:ix. lik.|(^-)/^+3/2' 

where B[dx) is the Fourier transform of a Brownian measure, 



^^^^ ~y_oo |j;|i+2(^^-i)/^ 
and the multiple stochastic integral is in the sense of [8]. 

For H G (1/2, 1) and S N* given, the Hermite process defined by (4.3) 
was studied independently in [9] and [25]. Its increments are stationary and 
its covariance is 

{ZI)W^{Z2)] = i(klP^ + |Z2P^ - \ZI - Z2P^). 

It is self-similar and i?-Holder. It is Gaussian if and only if A' = 1; thus, it 
is a fractional Brownian motion if and only if = 1. As a consequence, the 
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result of [19] corresponds to the case of = 1 in Theorem 4.1. Moreover, this 
result is in dramatic contrast to the short-range case where the asymptotics 
does not depend qualitatively on 

Proof of Theorem 4.1. Following [9] or [25], we find that the finite- 
dimensional distributions of the antiderivative of converge to those of 
W^; therefore, Ai is satisfied. Next we show that A2(2) and A3(2) hold. 
Because of the stationarity of m it is enough to show that 

(4.4) ¥.[v{Q)v{z)]^CuZ-^'' asz^oo 

for some constant > 0. By the Hermite expansion we can write 



oo 



Using the properties of the Hermite polynomials we get 

k=K 



(4.5) 



k=K 

Therefore, we need to study the limit of 



oo 



.^^E[Ko)K^)] = E 

k=K "•'^O 

Observe that for k = K we have z'^^rm{z) ~ c as z — t- oo, and for k > K 
we have z'^^rm{z)^ — >■ 0. Moreover, we have the uniform upper bound for z 
sufficiently large 

'^^^^%^^|r™(z)|'=<"^^^^' 



Using the fact that X^^^^ '^^W' ^ (^•^) follows from the uniform conver- 
gence theorem. □ 



5. Application to multifractal media. In this section we study the case 
where the asymptotic medium is described in terms of a multifractional pro- 
cess. In all the situations described above, the media were asymptotically 
expressed in terms of fractional processes. A drawback of fractional processes 
for applications is the strong homogeneity of their properties, which are de- 
scribed by their (constant) Hurst index. Therefore, multifractional processes 
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have attracted much attention [3, 21]. Multifractional processes have locally 
the same properties as fractional processes. Their properties are governed 
by a (0, l)-valued function h which is called the multifractional function. 
Some of the main properties are that multifractional processes are locally 
self-similar, and their pointwise Holder exponents vary along their trajec- 
tory. In particular, multifractional processes are relevant in order to describe 
nonhomogeneous media. Before stating the main result of this section, we 
mention that the most famous multifractional process is the multifractional 
Brownian motion. It was independently introduced in [3, 21] and can be de- 
fined from the harmonizable representation of fractional Brownian motion 
for every z 

(5.1) = _|^£_Z_B(,,,, 

where B is the Fourier transform of a real Gaussian measure B, and the 
constant C{H) is a renormalization constant and can be written as 

/oo \ -ix _ 1 |2 _ 

Now we consider a (0, l)-valued function /i, and we substitute H by h{z) for 
every z to obtain 

where the constant C{z) is a renormalization function. 

We shall here use a different framework for the multifractal modeling that 
is convenient for the asymptotic analysis and describe this next. We assume 
that z/^ has the form 

u^{z)=e<^^-^u(^^,^ forzG[0,Z], 

where k is a positive function, and is a field that is written as 1^(2:1, 22) = 
^{m{zi^h{z2))) for every zi and Z2 where: 

• <I> is a continuous function with Hermite index 1. 

• /i is a continuous function taking values in [/i_,/i+] C (1/2, 1). 

• m = {m{z,H)}z,H is a centered and continuous Gaussian field such that 
E[m(2;, -ff)^] = 1 for every z and H and such that there exists a continuous 
function r: [/i„,/i_|_]^ — )• (0,oo) (that we call the asymptotic covariance of 
m) such that 

lim sup \r[Hi,H2) 

^1-^2^°° (Hi, /fa) 

- {zi - Z2f-"'-"'nm{zi,Hi)m{z2.H2)]\=Q. 
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These assumptions describe that the field m has the long-range property 
with respect to the variable z. They also express that for each H, the process 
m{-,H) is stationary and asymptotically fractional because it satisfies the 
classical invariance principle. As established in [6] this field enables us to 
define a process that is asymptotically multifractional. 
Applying Theorem 3.1 we now get the following theorem. 

Theorem 5.1. Let 7(2) := r - k{z), and assume h{z) = (2 — j{z))/2. 
Then, as e goes to 0, {a'^(Z, s)}^ converges in distribution in the space of 
continuous functions endowed with the uniform topology to the random pro- 
cess {a{Z,s)}s that can he written as 

(5.4) a{Z,s) = f{s-\Sh{Z)), 

where Sh is a centered Gaussian process with covariance for zi , 22 > given 
by 

(5.5) E[Shizi)Sh{z2)] = J{lf dui du2n{uuU2), 

Jo Jo 

where 

n{ui,U2)= n{ui ,U2-MUI)MU2))\UI-U2\ M«i)+M-2)-2 

with 

(5.6) nizi,Z2;Hi,H2) = r{Hi, H2)U,>z, + v{H2,Hi)l,,<,,. 

The process Sh was introduced in [6]. This process is continuous and 
multifractional in the sense that its pointwize Holder exponent is /i(to) at 
the point to: 

sup< H, hm —rjj = > = h{tQ). 

Notice that in the case of h is constant Theorem 5.1 corresponds to the 
result of [19]. 

Proof of Theorem 5.1. By the same procedure as in proving (4.4), 
we get from the asymptotic assumptions for {m(z,H)} that 

hm sup \{zi-Z2f-'''-''mHzi,Hi)uiz2,H2)] 

-J{lfr{Hi,H2)\=G. 

If we denote, respectively, v^ and w'^ the antiderivatives of z 1— >• v^{z) and 
z I— )■ e^'*^^)~^m(z/e^, /i(z)), then, by using the same argument as above we 
also get 

(5.7) limE[|?;^(z) - J {l)w^ {z)\'^] = 0, 

£— >0 
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which imphes that the convergence of the finite-dimensional distributions of 
can be reduced to those of w^. Hence, without loss of generality and from 
the point of view of the analysis we can assume that $ = Id and work with 



u%z) = J{l)e 



,h{z] 



Following [6] , the finite-dimensional distributions of the antiderivative of z/^ 
converges to those of S^, and thus Ai is satisfied. Now we check A2(2). We 
let 5 > and thanks to the asymptotic assumption on m, there exists 
such that for every zi, Z2 and e satisfying |zi — Z2I > £^-^5 we have 

2-H1-H2 



sup 

[Hi,H2) 



Zl - Z2 



E 



Zl 



n 



Z2 



^1 ^2.^ rr 



<5. 



Then, noting that Tl{zi/e'^ , Z2/£'^ , Hi, H2) =Tlizi, Z2, Hi, H2) and substitut- 
ing {Hi,H2) by {h{zi),h{z2)) we get 

2-h{zi)-h{z2) 



Zl - Z2 



E 



Zl 



Z2 



m\ -^,h{zi) ]m\ -^,h{z2) 



n{zi,Z2;h{zi),h{z2)) 



<5. 



Letting TZ*{zi,Z2) :=TZ{zi,Z2;h{zi),h{z2)) and noticing that sup(l/7^*) < 
CO (because inf IZ* > 0) we obtain 

\E[u%ZiX{z2)]-n*{zi,Z2)\zi-Z2\''^'^^+'''^''^-^\ 
< Sn* (zi ,Z2)\ZI- ^2 |M-i)+M^2)-2 sup(l/7^* ) ^ 

which proves A2(2). It remains to check A3 (2). Let /3 > 0. Because of the 
boundedness assumption on m, there exists a constant Ci{p) > so that for 
every zi, Z2 and e satisfying \zi — Z2\/ < p, we have 

' Z2 



E 



m 



Zl 



Hi ]m 



<Ci{p). 



Thus, 



\E[v%zi)v'{z2)]\ < C7i(p)e2M-i)+2M^2)-4 

= Ci{p)\zi - Z2\''^''^+''^'^^-^ 



Zl - Z2 



2-h{zi)-h{z2) 



< Ci{p)\zi - 22|'^(^l)+'^(^2)-2p2-/i(^i)-/i(22) 
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where C2{p) can be chosen such that Ci(p)p2-'^(^i)-'^(^2) < C2{p). So A3(2) 
is satisfied and the proof can be concluded by applying Theorem 3.1. □ 



We finish this subsection by applying Theorem 5.1 to an example that 
was mentioned in [6]. Let us consider Wh defined as in (5.1). We let 



(5.8) 



miz,H) = WHiz + 1) - Wniz). 



We compute the covariance between m{zi,Hi) and m{z2,H2) for every zi, 
Z2, Hi and H2 



(5.9) 



E[m{zi,Hi)m{z2,H2)] 
^l C{{Hi + H2)/2f 
2 C{Hi)C{H2) 



\Zl - Z2 



IH1+H2 



1 + 



1 



Zl - Z2 



H1+H2 



+ 



Zl - Z2 



H1+H2 



By Taylor's formula we get that the asymptotic covariance r of {m(z, H)}z,H 
can be written as 

(5.10) r{H,,H2) = \{H,+H2){H,+H2 - l)^^|^^y^f^- 

Then applying Theorem 5.1 we get that {a'^{Z, s)}s converges in distribution 
to a{Z, s) = f{s — ^Sh{Z)) where 



(5.11) ShiZ) = J{l) 



-ixe 



C(/i(u))|2;|M«)+V 



-^du\B{dx). 



As mentioned in Section 6.1 of [6], we also can observe that if we assume 
that h is differentiable then we can write Sh{Z) as 



-iZx 



Sh{Z) = J(l) / B{dx 
(5.12) 



1 



lC(/i(Z))|x|'^(^)+i/2 



|x|M«)+l/2 \C{h{u)) 



C'{h{u)) \^,, 



which means that Sh{Z) is the sum of a multifractional Brownian motion 
as in (5.3) and of a regular process. 
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6. A non-Gaussian and multifractal medium. In this section we study 
the case of a medium that generalizes the media discussed above. We define 
{m{z, H)}z^H for every z > by 

(6.1) m{z,H) = —^ I exp{-izx)ip{x)\x\^/'^-^ 13{dx), 

where H € (1/2,1), C{H) is a renormalization constant, V is a complex- 
valued symmetric function and B{dx) is the Fourier transform of a real 
Gaussian measure. We assume that ij: is continuous, "0(0) = 1 and satisfies 
|'(/'(x)| = 0\x\^Qo{\x\~^)- Notice that the family of processes defined by (5.8) 
in terms of fractional Brownian motion {Wh{z)} z,H is an example of such 
a process. 

Thus, {m(z, H)}z^H is a centered Gaussian field and its covariance can be 
written as 

/ N / M f e^V>(i(z2 — zi)x)\^(x)\'^ , 
(0.2) nm(.„H,)mU,.H,)] = 

Now we consider a function h that takes its values in [/i„, ^,4.] C (1/2, 1) and 
a truncation function <I> with Hermite index K gN* . We define as 

i.-(2)=e«W-z.(^^,2 
where 

^{zi,Z2) = ^{m{zi,hK{z2))) 

with 

7 , N Kz) - 1 , 
hK{z)= +1. 

We can then show that satisfies assumptions A2(2) and A3 (2). In partic- 
ular, we have 

(6.3) ^\V^{Z{)V'{Z2)\ ~ '^-^^v(hK{z{)~hK{z2))W " Z^t'^'^^^'^^'^^-'' 

when \z\ — Z2\/s'^ goes to 00 assuming that k,{z) — t = 2h{z) — 2, and r is 
defined as in (5.10). Therefore, because Theorem 3.1 says that, under long- 
range assumptions, the asymptotic behavior of a^{Z,s) is essentially given 
by the limit of v^{z), we can conclude by the following result. 

Theorem 6.1. As e goes to 0, {a^{Z,s)}s converges in distribution in 
the space of continuous functions endowed with the uniform topology to the 
random process {a{Z,s)}s that can be written as 

(6.4) S(Z,s) = /(s-i5f(Z)), 
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(6.5) 

where 



f ^ 



gh,K{z,Xi,...,XK) 



-IXk 



j(K)e-™ELi^fc 
K\C(hK{u))^ |xfc|''^(")+V2 



K 

n 



■ du. 



Notice that the process is equal (in distribution) to Wj^ of Section 
4 if /i is a constant equal to and is equal to Sh of Section 5 if = 1. 
Because of these facts, is in general non-Gaussian and multifractional. 
This shows that under general long-range assumptions the asymptotic time- 
shift is neither Gaussian, nor homogeneous. This is in dramatic contrast to 
the short-range case where the time shift is a Brownian motion, which is 
homogeneous and Gaussian. 



Proof of Theorem 6.1. We let 



v'{z) 



z/^(u) du ■ 



m 



and 



wk{z) 



due''^^-^-'PK[m{^,hK{i 



Using the same arguments as for the beginning of the proof of Theorem 5.1, 
and the fact that the Hermite index of ^ is K, we get 

2-\ 



(6.6) 



limE 

£->0 



JjK) 
K\ 



wk{z) 



0. 



Then using the formula (see [14], for instance) 



Pk[ / ^{x)B{dx) 



K 



4){xk)B{dxk) 



k=l 



for every G -L^ 



we get 

g2/i(n)-2 



K 



du- / e-^^i^i"^^/^ TT 

Cih{u))K J^K J-^i \^^\h^{u)-l/2 



B{dxk) 



K 



du W 



ipixk) 



2hiu)-2 ^ 

B[dxk) l ^^ e-^"^.^=i^^/" 



C{h{u)) 



K 
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Then we make the substitution — )• e'^Xk for every k 
We let 

= f fduf] -^t^^^B{dxu)—^e-^^^U^K 
LkJo |rEfc|'^^(«)-i/2 ^ '''C{h{u))K 

The self-similarity of the Brownian motion gives that B[e'^ dxk) is equal in 
distribution to eB{dxk), then we get that 

e f.d.d. ~s 

where ^'='^' means the equality of the finite-dimensional distributions. Then, 
using the assumptions on ip, we obtain the convergence a.s. of the finite- 
dimensional margins of ^^w^k those of , and thus the convergence of 
the finite-dimensional distributions of to those of , so Ai is satisfied. 
Now, as observed at the beginning of this section, using (6.2) and by the 
same procedure as in the proof of Theorem 5.1 we show that A2(2) and 
A3(2) hold. We can then conclude by Theorem 3.1. □ 



7. Numerical illustration. We illustrate our results with some numerical 
simulations. In order to show the differences between the mixing and the 
long-range cases, numerical simulations of the transmitted pulses centered 
around the travel time are presented in [19]. They are carried out with a 
fractional white noise medium with Hurst index if = 0.5 (corresponding to 
the mixing case) and H = 0.6 (corresponding to the long-range case). These 
examples illustrate that the pulse shape is not affected by the random fiuctu- 
ation of the medium when H = 0.6 and that it is modified via a convolution 
with a Gaussian kernel when i/ = 0.5. 

Here we aim to illustrate the differences between fractional and multifrac- 
tional cases. We present simulations of the asymptotic travel times we obtain 
for media with long-range correlation and different multifractional functions. 
For the sake of simplicity we restrict ourself to the Gaussian case presented 
in Section 5, and we let the propagation distance be one. For a fixed multi- 
fractional function /i, the method we use to simulate the asymptotic travel 
time Sh is based on the method presented in [22] (pages 370-371) and the 
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invariance principle proved in [6]. We first simulate the fractional white 
noise {Yj{H)}j of index He (1/2,1) as in equation (7.11.1) of [22] (page 
371). Then, using Theorem 2 of [6] we can use Yl^^li N-'''^^/^^Yj{h{j /N)) 
to approximate Sh{t). In Figure 1 we show a trajectory of Sh with an in- 
creasing multifractional function. In Figure 2 we show a trajectory of Sh 
with an periodic multifractional function. In both figures we can observe 
that the regularity varies along the trajectory according to the local Hurst 
index. Modeling of this kind may, for instance, be relevant in the case the 
multiscale crust of the sedimentary earth or in the context of the turbulent 
atmosphere. In both cases the field is typically strongly anisotropic with a 
roughness that depends on depth or height, respectively. 



8. Proof of Theorem 3.1. We first give an outline of the proof. As recalled 
in Section 2 the process {a^{Z, s)}s can be written in terms of the propagator 
P^, and thus the study of the convergence of {a^{Z,s)}s can be analyzed 
via asymptotic properties of P^. The propagator P^ satisfies the equation 



Multifractional function h Trajectory of S_h 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Fig. 1. Trajectory of Sh with an increasing multifractional function. 
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Multifractional function h 



Trajectory of S_h 




I I I I I I I I I I I I I I I I I 
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 




I I I I I I I I I I I I I 
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



I I I 
0.8 0.9 1.0 



Fig. 2. Trajectory of Sh with aperiodic multifractional function. 



that we can write in the form 



(8.1) 
where 

and 



1 
-1 



-1 

1 



i 

1 



and vf, V2 and are three processes of bounded variation that we can write 
as 

vl{z)= [\%z')dz', 
Jo 

vl{z) = u''{z')cos(^uj^\dz'. 



v%z')sn\(^uj^ dz. 
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Thanks to T. Lyons's rough paths theory for which we recah some tools in 
the Appendix we shall see that the convergence of can be reduced for a 
convenient topology to the convergence of the process v*" defined as 

yr^ ■.= [v\,vl,vl). 

Hence, we first prove the convergence of v^, then by Theorem A.l (see the 
Appendix) we deduce the convergence of in Section 8.1 and thanks to 
(2.13), we finally conclude by the convergence of {a^{Z,s)}s in Section 8.2. 

8.1. Convergence of the propagator. Using Theorem A.l and the expres- 
sion (8.1), the asymptotic study of the propagator is reduced to finding 
the limit in a rough path space of v"^ := {vi^v^^v?^). This is the aim of the 
following lemma. 

Lemma 8.1. There exists 7* G (0, 1) such that for every p > 2/(2 — 7*), 
as E goes to 0, the increments of converge in Vtp to those of V which can 
he written as 

v = (y,o,o). 

The proof of Lemma 8.1 is based on establishing several technical lemmas 
that we do next. We let r,^e(x,y) = E[i/^(x)z/^(y)]. 

Lemma 8.2. There exist C and 7* so that 

\rye{x,y)\ < C|x 

for every x and y. 

Proof. The assumptions of Theorem 3.1 imply that for every S > 
there exists zs > such that for — 7/| > e^zs we have 

(1 - 6)R{x,y)\x - y\-^(^'y) < r,s{x,y) < {1 + 6)R{x,y)\x - 

Hence, taking 5 = 1 we get that for \x — y\ > e^z\ we have 

< riye(x,y) < C\x — 

Moreover, thanks to the assumptions of Theorem 3.1, we know that there 
exist C^j and 7^^ so that for — y| < e"^zi we have 

0< |r^=(2;,y)| < C^Jx - . 

By choosing 7* := max(7+,72;^) we get that there exists 7* so that 

\E[u'{x)u'{y)]\<C\x-y\-^* 

for every x and y. □ 



20 R. MARTY AND K. S0LNA 

Lemma 8.3. For every z G [0, Z], as e goes to the sequences V2{z) and 
v^{z) converge to 0. 

Proof. Without loss of generality we present the proof only for ^2(2) 
and with 2uj = 1. We have 

^[v2izf]= j dx j dy cos (^-^^ cos(^^^r^e{x,y) 
= Il{z)+II{z), 



with 



niz) = dx dy cos cos Rix, y)\x - , 

Iliz) = J^'dxJ^'dy cos(^^) ^°^(|^) (^-^(^'y) - " yl"^^"'^^)- 

Let 6 > 0, and because of the assumptions of Theorem 3.1, we have that for 
\x — y\ > e^Zf, (with 25 sufficiently large) |rjye(x,y) — R{x^y)\x — y\~^^^^y)\ < 
6R{x, y)\x — y|~'>'(^'f) for every e. Combining this with Lemma 8.2 we obtain 

\II{z)\<6 Tdx r dyR{x,y)\x-y\-^^-'y^ 
Jo Jo 



Cs dx dy\x — y\'^l 
Jo Jo 



so that 

limsup|/|(2:)| < dx dy\x - y\~'^^^'y\ 

e->0 JO Jo 

The inequality above is valid for every 6 > 0, and we conclude 

lim/|(z)=0. 

£-5>0 

We can deal with Ifiz) using a Riemann-type result. Indeed, the function 
R: {x,y) I—)- R{x,y)\x — y\~'^(^'y^ is integrable on = [0,2;]^, so we can ap- 
proximate it by a sequence of constant by step functions {Rn)n such that 

lim dx dy\R{x,y) - RNix,y)\=0. 

N^OoJq Jq 

Moreover, we can write 



dx j dy cos (^-^^ cos (^■^^RN{x,y) 



+ dx dy\R{x,y) - RN{x,y)\ 
Jo Jo 
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for every e and N. We easily see that 



hm / dx dy cos 
e^oJo Jo Ve 

so that 



Umsup|/f (z)| < dx dy\R{x,y) - RN{x,y)\ 
e->0 JO Jo 

for every A^. This finally shows 



and then 



lim If (z) = 



limE[t>|(z)^]=0, 



which completes the proof. □ 

Now we deal with a technical lemma regarding the increments of v^. 

Lemma 8.4. There exist C > and 7* £ (0, 1) such that for every z, C 
and e > we have 



E[K(z)-v^(C)f]<C|z-C| 



2-7* 



Proof. Because of Lemma 8.2 there exists 7* so that |E[z^^(x)z/^(y)]| < 
C\x — y\~'^* for every x and y. Then, for every j = 1,2,3, we have (taking 
z>C) 

n\v'jiz)-v'j{C)\^]< J\x J^' dy\E[u'{x)i^%y)]\ 



<Cy dx J dy\x — y\'^ 

-(i-r)(2-r)' ^' ' 

which completes the proof. □ 

In the sequel we shall use the notation i?^, := (2 — 7*)/2. Using the above 
lemmas we next deduce the following lemma which deals with identification 
of the limit. 



Lemma 8.5. The process V defined in Lemma 8.1 is a.s. continuous (up 
to a modification). Moreover, as e goes to 0, v"^ converges to V in the space 
of continuous functions endowed with the uniform norm. 
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Proof. Assumptions and Lemma 8.3 give the convergence of finite- 
dimensional distributions of to those of V. Using then the Kolmogorov 
criterion [4], Lemma 8.4, and the fact that 2H^ > 1 we get the tightness of 
(v'^)^ in the space of continuous functions endowed with the uniform norm 
which estabhshes the proof. □ 

Thanks to Lemma 8.5 we conclude with the proof of Lemma 8.1 by es- 
tablishing the tightness in a rough paths sense. 

Lemma 8.6. For every p > l/H^, the sequence (v^)^ is tight in Qp and 
the process V is a.s. of finite p-variation. 

Proof of Lemmas 8.1 and 8.6. Let q G {1/H^,p).ln view of Lemmas 
A.l and 8.5 it is enough to prove 

(8.2) lim supP[yg(v^) >A]=0. 

A->+oo £>o 

Using Chebyshev's inequality, the fact that q <2, Lemma A. 2, the Holder 
inequality and Lemma 8.5 we find 



n=l ^ ^ 

and since qH^ > 1 we deduce (8.2). □ 



A^ 

„ +00 

n=l k=l 
^ +00 2" 

n=l k=\ 

n=l k=l ^ 



Finally, we can now derive the following lemma which deals with the 
convergence of the propagator. 

Lemma 8.7. Let {tJi, . . . ,t<;„} to he a collection of frequencies. Then, as 
e goes to 0, the propagator vector (P^^ , • • • , ) converges in distribution in 
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the space of continuous functions to {P^^j^, . . . , P^^^) which is the asymptotic 
propagator P^^ that we can write as 



Proof. By combining Theorem A.l, (8.1) and Lemma 8.1 we get that, 
as £ goes to 0, converges in distribution in the space of continuous func- 
tions (endowed with the topology of the uniform convergence) to the solution 
Poj of the following system of equations: 



This concludes the proof. □ 

We remark that the situation here contrasts with the short-range case. 
Indeed, the asymptotic propagator is driven by one process in the long-range 
case whereas it is driven by three processes in the short-range case. 

8.2. Conclusion of the proof The remaining part of the proof of Theorem 
3.1 follows the lines of [5, 10]; however, we present it here for completeness. 
Recall that thanks to the formula (2.13) we can write a^(Z,s) in a Fourier- 
type formula using the transmission coefficient 



with the transmission coefficient being a functional of the propagator P^. 
We shall use Lemma 8.7 to deduce the convergence of the transmitted wave. 
Let n G N, si < • • • < s.„ G [0,oo). We can write: 





(8.3) 




E[a%Z,si)---a%Z,Sn)] 




n 




xE[T^^{Z)---T^^{Z)]du;,---dio, 



Thanks to Lemma 8.7 we have that as e — ?• 
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E[a''{Z,si] 
1 



(2vr)^ 
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X E 



exp 



iV{Z) 



1 • • 



(27r)" J-J- 



nM-n^)/2) 

Li=i 

The tightness proof is similar to the proof of Lemma 3.2 in [5] and the 
convergence of a^{Z,s) follows. 

APPENDIX: DIFFERENTIAL EQUATIONS AND ROUGH PATHS 

In this appendix we fix p G [1, 2) and consider a closed interval / = [0, Z]. 
We define the variation of a continuous function w.I^ M" by 



fe-i 



i/p 



Vp{w):= ( sup^||w(zj+i) 



j=0 



where sup^j runs over all finite partition {0 = zq, . . . , Zk = Z} of / and where 
here and below || • || refers to the norm. The space of all continuous 
functions of bounded variation (1-variation) is endowed with the variation 
distance 

ll^llp = ^("w^) + sup 

ze[o,z] 

and is denoted by il.^ . The closure of this metric space is called the space of 
all geometric rough paths and is denoted by ilp. One of the most important 
theorems of rough paths theory is the following: 

Theorem A.l (T. Lyons's Continuity Theorem). Let^ G:M x M'' 



£(M,M'^) andF: 



be two smooth functions. Let y be the 



unique solution of the differential equation 

dy{z) = G{z, y{z)) dz + F{z, y{z))dw{z), y{z = 0) = yo, 



^Here £(R,R'*) [resp., £(R",R'*)] denotes the space of all linear maps from R (resp., 
') to R'*. 
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where w is a bounded variation function. Then ltd 's map X:w ^ y is con- 
tinuous with respect to the p-variation distance from r2^(M") to Vl'^{W^). 
Therefore there exists a unique extension of this map (that we still denote 
by I) to the space Qp{M"-) 

This theorem has been proved by T. Lyons and extensively studied and 
appHed (see [7, 15-17]). 

The proof of Theorem 3.1 is based on analysis of the tightness in the 
space of geometric rough paths. In the context of this we need to compute 
the p-variation for p> 1. To this effect we will need the following lemmas of 
which the first can be found, for instance, in [16], and the second in [15, 16]. 

Lemma A.l. Let q € [1,2) and (w^)e>o o- family of continuous random 
processes of finite q-variation whose associated family of probability measures 
is tight in the space of continuous functions on I and satisfying 

(A.l) lim supF[Vq{v') > ^] = 0. 

Then the family of probability measures associated to {v'^)e>o is tight in Jlp 
for every p> q. 

Lemma A. 2. For every n G N and every k = 0,1, . . . , 2", we let := 
Zk/2^. Let q G [1,2) and v be a function of finite q-variation. Then there 
exist two positive constants Ci, C2 which do not depend on v such that 

+00 2" 
n=l fc=l 
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